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Abstract 

Selection effects, connected with stochastic errors in source flux and 
threshold value determination are analyzed. Normal and normal logarith- 
mic distributions of stochastic deviations are considered. These two kind 
of distributions produce different effects on the source statistics. Appli- 
cations to Gamma Ray Burst statistics are discussed. A physical test for 
checking a close neutron star model of GRBs is suggested. 

1 Introduction 

Statistical investigation of samples of sources is a powerful method for analyz- 
ing their location, origin and evolution. Most important rcsiilts were obtained 
for distant radio sources (Longair, 1966). In combination with redshift mea- 
surements these data permit also to get estimates for cosmological parameters 
(Zeldovich and Novikov, 1975). 

The BATSE (Meegan et al,1992; Kouveliotou, 1994; Fishman and Meegan, 
1995) curve [log A'' — log(C/CTOi„)] gives very important information for making 
constraints on GRB models and understanding their nature. Nevertheless, it 
suffers from different selection effects, so it seems premature to use it for a 
critical choice of GRB models. This curve differs from the straight line with a 
slope 3/2, and the observed isotropy of GRB distribution on the sky is consistent 
with the following models: 

1) Nearby neutron stars from the disc population. 

2) Galactic halo neutron stars. 

3) Cosmological model with bursts coming from sources with rcdshifts z ~ 
1^2. 

The last two models also explain deviations from the 3/2 line. The third one has 
extreme demands (neutron star collision); huge energy release in soft gamma 
region with no counterpart in optics or radio. While some theoretical models 
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based on fireball expansion (Mcszaros & Rees,1993) seems available to repro- 
duce the main properties of GRB, it is questionable to get such a fireball in 
neutron star collision (Isern et.al.,1995). There is also a problem of extended 
hard gamma - emission, accompanying the main burst (Hurley et.al.,1994). The 
halo model is posing restrictions (Hakkila et al,1995; Bulik and Lamb, 1995) to 
the properties of the neutron star kinematics (speed at the origin ~ 10'^ km/s), 
GRB fluence (narrow strip around 10'*^ ergs) and location (at the outer edge of 
a sphere with R ~ 350 kpc). 

Here we investigate the influence of stochastic errors on the shape of the 
[logiV — \og{C /Cmin)] curve. Such type of errors have been first taken into 
account by Eddington (1913,1940), who considered small stochastic errors in 
determination of the observed magnitude of stars. Fluence or peak flux of 
GRBs are determined with much larger errors, related to the uncertainty in 
the burst angular direction, its spectral and temporal variability, background 
fluctuations. 



2 Source statistics in presence of stochastic er- 
rors 

Let X be a number density of events, registered by an observer with a flux C . If 
V IS a, constant burst frequency per unit volume, L is a constant peak luminosity, 
Cx is a variable threshold flux, g{Cx) is a threshold distribution function, then, 
following Schmidt et al (1988), Band (1992), Petrosian (1993) we may write 

X{C) - Airy j 5{C-^)e{C-Cx)g{CxydrdCx (1) 

a) Logarithmic normal distribution. Consider first a classical case 
where threshold effects are neglected. Actually, Eddington (1913,1940) con- 
sidered a case with no threshold, when all sources could be registered. For that 
case we have instead of (|l|) 

^{^O ^ A^v j 5{C-^)r^dr (2) 



which gives 



x(C)^f:&V"|S;, (3) 



3/2 

An) C5/2 

Following Eddington (1913), take errors, according to a normal Gaussian distri- 
bution over log C. The observed flux C could be produced by the source, whose 
real flux is C" with the probability 



-r-^exp 



(logC"-logC)2 
A2 



d{logC') (4) 
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The observed number density of events x(C) is connected to the real distribution 
X(C') by the relation 



xic) 



1 



xx(C") exp 



Using (H) in (||) we get after integration 



A0F 

(logC"-logC)2 



(5) 



rf(logC') 



X{C) ^x{C)eT^^ 



(6) 



So, Gaussian logarithmic statistical errors in absence of threshold does not 
change the slope of the number density curve, increasing it by a constant coef- 
ficient (Eddington, 1913) 

b) Normal logarithmic distribution with variable threshold. For 

a variable threshold with g{Cx) ^ <5(Co) we need to consider a distribution 



(Schmidt et.al. 



over ^ 



c 



(7) 



x0(eC,-C,)5(a)r2 drdC, 
which gives after integration 
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Consider now a Gaussian logarithmic distribution of errors like with ^ and 
^' instead of C and C . While it is impossible to register events under the 
threshold, we shall use the interval 1 < ^' < oo for possible real values. Then 
we have instead of (H) 



m = ^ 



xiC) cxp 



(log^'-log^^ 



d(log^') 



(9) 



with 



B = 



cxp 



{x - logQ^ 

A2 



dx 



Using (H) in (0) we get 
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Here 

/■oo „ /■(log5)/A 

B^A{ e-'- dz + e-'- dz) (11) 

^0 "'0 

where Erf{x) — e^^^ dz. An observed distribution S(^) obtained by 

integration of (hw over is written as follows 



oo 

(12) 

4tt J x^l^ 



With 



The function logS as a function of logf is represented in fig.l for A = 1. The 
corresponding curve for A = 0.1 is cannot be distinguished from the straight 
line with slope 3/2 indicated there. Contrary to the case with no threshold 
(Eddington,1913), where stochastic errors do not change the slope of the curve 
[log TV — log C] , stochastic errors in presence of threshold may considerably de- 
crease the slope of the curve in the vicinity of the threshold. It happens, because 
the threshold is acting like a border, which cannot be crossed by the bursts from 
both sides. So faint bursts appear more like stronger ones, and the spreading of 
stronger bursts is almost equal in both directions. 

c) Normal distribution. While we do not know exactly what law is de- 
termining stochastic errors, let us consider also a normal distribution of errors 
around the value ^ itself. While our eye has a logarithmic response to the signal, 
some X - ray counters are proportional. May be the true distribution is even 
more complicated and does not follow a Gaussian law neither in logarithms nor 
in the values themselves. 
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Figure 1: The curve [log iV— log C {max) / C{thr)\ in presence of stochastic errors, 
distributed according to normal logarithmic distribution; 1 - straight line with 
a slope 3/2, corresponding to A = 0; 2 - curve with A = 1; C{max) is the peak 
intensity of the burst; C{thr) is a corresponding threshold value. 
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Assume that the burst, registered with the ratio ^ — C/Cx may in reahty 
correspond to the ratio ^ with the probabiHty 



exp 



(14) 



Then taking into account only the events over the threshold we get instead of 
(i 



with 



poo 


)exp 




poo 










A? J 



^2^ 



Taking into account the (pf) we get 



where the integrals may be expressed as 



L \ 2t:vA h 



/i = / y ^/^expU [y~^] \ dy 



Ai 



Ai 



For the distribution we are looking for we get 



3/2 



with 



Iijx) 



Di{x)dx 



(15) 
(16) 

(17) 

(18) 
(19) 

(20) 
(21) 



For sufficiently large values of ^, when the function Erf{x) is very close to unity, 
the expression for S(^) in (EQ) may be written as 



3/2 



3/2 



X e ^ dz + 



Ai 



Ai 



-3/2 



(22) 
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The plot of log ^(^) as a function of log^ is represented in fig. 2 for Ai = 1, 10, 
together with the straight line with the slope 3/2. 

3 Discussion 

BATSE consists (Fishman, 1992) of eight detectors, arranged on corners of the 
Compton Gamma Ray Observatory (CGRO). Bmst registration is done by large 
area detectors (LAD), optimized sensitivity and directional response. The eight 
panels of LAD are parallel to the eight faces of a regular octahedron. Since 
a regular octahedron is comprised of four sets of parallel intersecting planes, 
every detected burst will be viewed by four detectors. LAD are sensitive in 
the energy range 20-600 keV. The burst is registered, when 5.5a excess over 
17 s background rate is registered at least by two detectors (Fishman et al, 
1992). The background in LAD in the burst trigger energy range 60-300 keV 
varies between approximately 1500 counts/s and 3000 counts/s per detector 
during most portion of the orbit and above geographic latitudes of about 22° 
the background increases considerably. 

Each event with the ratio ^ of the peak luminosity to a local background is 
detected with an error due to the following circumstances. 

1) Spectral dispersion of GRBs cannot guarantee that a peak value in the 
BATSE spectral region (BSR) is equal to a real one. It may be several times 
larger if the region of maximum radiation lay outside BSR. 

2) Angular dependence of the detector sensitivity, especially in the region of 
a steep dependence around 50° (Fishman, 1992) imply errors in determination 
of ^ because of poor angular localization of the source. 

3) Different time duration of bursts lead to nonuniformity of the source 
sample, where a same peak luminosity may be related to bursts with a total 
flux or average luminosity varying by orders of magnitude. Conversely very 
different peak luminosities may correspond to bursts with the same total flux. 
This would imitate stochastic errors of the same order. 

4) Nonuniformity of GRB detection conditions, when an event may be reg- 
istered by 2 or 3 or 4 detectors determines an additional source of dispersion. 
So, uncertainty in the ^ value represented by the dispersion equal to 10 threshold 
levels, described above, does not seems to be overestimated. 

For Al = 10, corresponding to A = 1 in the case of logarithmic normal 
distribution of errors, the shape of the curve is changed considerably. The 
changes become noticeable at ^ ^ 30 (see fig. 2), which is much larger, than the 
average dispersion value Ai = 10. In the case of logarithmic dispersion the 
influence of stochastic errors starts at ^ approximately equal to the value of the 
average dispersion (see fig.l). Results, represented in fig. 1,2 illustrate the large 
importance of different stochastic errors in the form of the [log N — log{C / Cmin)] 
curve and are not intended to explain directly the corresponding BATSE curve. 

We have not taken into account other kinds of threshold influence, leading to 



8 



Figure 2: Same as in Fig.l for normal distribution; 1 - straight line with slope 
3/2, corresponding to Ai = 0; 2 - curve with Ai = 1; 3 - curve with Ai = 10. 
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additional deviations from the 3/2 slope (Hartman and The, 1993; Lingenfelter, 
1995). Combination of these effects must be taken into account in analyzing the 
[log TV - \og{C /Cmin)] curve for the BATSE data sample. 

Let us note, that stochastic errors from the photon noise, to which Eddington 
(1913) applied his calculations, are much smaller than possible errors in BATSE 
data caused by above mentioned reasons. Even in statistical analysis of G- stars, 
situated uniformly inside the Galactic disc, the value of < V/Vmax > which must 
be equal to 0.5 for uniform unbiased sample, falls down considerably around two 
threshold flux, which is connected most probably with a loss of faint sources 
(Harrison et al, 1995). 

In view of this situation it seems preliminary to wipe away a close Galactic 
origin for GRBs and all variety of models (see e.g. Ho et al, 1992) should remain 
under discussion. 



4 Testing the physical origin of GRBs 

If GRBs are connected with starquakes on nearby neutron stars (Bisnovatyi- 
Kogan et al, 1975) it is worth (Bisnovatyi-Kogan, 1993) to monitor close young 
pulsars (Geminga) for catching the moment of the quake and comparing it with 
GRB search data. If the hard tail of GRBs (Hurley et al, 1994) is connected with 
the excitation of submsec proper oscillations of the neutron star after starquake, 
and hard gamma ray emission in GRBs is produced by the same mechanism 
as in radiopulsars (Bisnovatyi-Kogan, 1995), one implicit test for GRB origin 
may be suggested. If starquakes in radiopulsars lead to the excitation of such 
oscillations, they can lead to the appearance of resonant modes in a frequency 
spectrum of radioemission. 

The electrical field, generated on the surface of an oscillating neutron star, 
is of the order of (Muslimov and Tsygan, 1986) 

^o.c-^B«^^B«^^B (23) 

c R c c 

When the amplitude of oscillations is large enough 

5R Qrot ^ 10-4 

this oscillating field could modulate a pair cascade birth, leading to the appear- 
ance of a coherent high frequency mode in the frequency spectrum of radioe- 
mission. 

In old nearby neutron stars (silent ones) we may expect slower rotation and 
lower magnetic field than in radiopulsars. These neutron stars could become 
pulsars (in hard gamma as well as in radio) only temporally, for about a few 
hours after the quake, and produce a GRB, if the electrical field, induced by 
oscillations is higher than the threshold field for pair cascade generation. 
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The best object for such testing is the Vela radiopulsar, where strong ghtches 
are observed almost every year (MacCuUoch et.al.,1987). It is necessary to be 
able to make a frequency analysis of the radio data very soon after the quake 
for checking the existence of resonance frequencies with periods less then one 
millisecond. If the model of close Galactic GRB with its logical consequence 
listed above is true, we may expect to sec high frequency resonance oscillations 
only during a limited period of time of the order of 90 minutes after the visible 
glitch. Radio observations of Vela soon after glitch with high time resolution 
and accurate frequency analysis could be more informative than gamma ray 
observations (Hartmann et al, 1992). 
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